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We analyze the equilibrium micromagnetic domain wall structures encountered in Permalloy strips 
of a wide range of thicknesses and widths, with strip widths up to several micrometers. By per¬ 
forming an extensive set of micromagnetic simulations, we show that the equilibrium phase diagram 
of the domain wall structures exhibits in addition to the previously found structures (symmetric 
and asymmetric transverse walls, vortex wall) also double vortex and triple vortex domain walls 
for large enough strip widths and thicknesses. Also several metastable domain wall structures are 
found for wide and/or thick strips. We discuss the details of the relaxation process from random 
magnetization initial states towards the stable domain wall structure, and show that our results are 
robust with respect to changes of e.g. the magnitude of the Gilbert damping constant and details 
of the initial conditions. 

PACS numbers: 75.60.Ch, 75.78.Cd 


I. INTRODUCTION 

During the last decade, a lot of effort has been de¬ 
voted to understand static and dynamic properties of 
magnetic domain walls (DWs) in ferromagnetic nanos¬ 
tructures such as nanowires and -strips. These studies 
have been largely driven by promising technological ap¬ 
plications based on domain walls and their dynamics, in 
particular memory 1 2 ^and logic devices 3 5 . In typical ex¬ 
periments DWs are driven by either applied magnetic 
field! 6 7 or spin-polarized electric currents 8 11 . The re¬ 
sulting DW dynamics depends crucially on the micro- 
magnetic DW structure, typically involving various in¬ 
ternal degrees of freedom. These are essential e.g. for 
the emergence of the Walker breakdowrP^, an instability 
occurring when the DW internal degrees of freedom get 
excited by a strong external drive (a magnetic field H 
or a spin-polarized current J exceeding the Walker field 
77w or current Jw, respectively), limiting the propaga¬ 
tion velocity of the DWs. 

Two main classes of ferromagnetic materials have been 
extensively studied within the strip geometry. Ma¬ 
terials with a high perpendicular magnetic anisotropy 
(PMAEEU) 

exhibit simple and narrow DWs of the 
Bloch and/or Neel type. For H > i7w or J > Jw, 
repeated transitions between these two structures are 
observed 16 ! The second class of system s includes soft 
(low anisotropy) magnetic materials 13 31 such as Permal¬ 
loy, where in-plane domain magnetization along the long 
axis of the strip is induced by shap e ani sotropy. By us¬ 
ing various experimental technique ! 17 * 18 ^ and micromag¬ 
netic simulations, it has been established that the equi¬ 
librium DW structures separating these in-plane domains 
are more complex, and depend crucially on the sample 
geometry®* 2 2 f Transverse DWs (TWs) and asymmetric 
transverse DWs (ATWs) are observed for narrow and 
thin strips®^, while in wider and thicker strips one 
encounters the vortex DW (VW)®- 24 ^ 26 ! 28 ^ . In addi¬ 
tion, various metastable DW structures with higher en¬ 


ergy may be found^L^EO. For H > 77w or J > Jw, the 
DW structures exhibit dynamical evolution: for TWs, 
repeated nucleation and propagation of an antivortex 
across the strip width takes placet Similarly, in VWs 
the vortex core performs oscillatory back and forth per¬ 
pendicular motion 19 . 

In Permalloy strips with even larger widths and/or 
thicknesses, one might expect also other, possibly more 
complicated equilibrium DW structures. For wider strips 
shape anisotropy is less important, implying that energy 
minima with more complex spin structures closing the 
flux more efficiently than TWs, ATWs or VWs may ap¬ 
pear. Indeed, e.g. double and triple vortex DWs have 
been observed in experiments on wide strips®, but they 
have been attributed to current-induced vortex nucle¬ 
ation resulting in metastable DW structures. Conse¬ 
quently, a pertinent and fundamental question is what 
are the possible intermediate equilibrium DW structures 
observable when the lateral Permalloy strip dimensions 
increase from those corresponding to the typical nanos¬ 
trip geometry (with TW, ATW or VW as the stable DW 
structure) to strip widths of micrometers and beyond. 

In this paper we present an extensive numerical study 
of the equilibrium and metastable micromagnetic DW 
structures in Permalloy strips, with the strip widths up to 
an order of magnitude larger than before 19 Contrary 
to previous studies focusing on comparing the energies 
of different a priori known DW structures®"®, we per¬ 
form micromagnetic simulations of relaxation dynamics 
from random initial states towards the stable DW struc¬ 
tures. In addition to the previously observed TW, ATW 
and VW DWs, we find also DWs with equilibrium double 
and triple vortex structures for wide and/or thick enough 
strips. The last structure is encountered only in the very 
largest system sizes we were able to simulate. Moreover, 
for wide strips we find a rich variety of metastable DWs 
with even more complex micromagnetic structures. We 
demonstrate that our results are robust with respect to 
changes of the magnitude of the Gilbert damping con- 
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FIG. 1. (color online) (a) Geometry of the Permalloy strip, 
(b) A top view of the magnetization in the initial state. Mag¬ 
netization points along the long axis of the strip within the 
two domains (as indicated by the arrows) forming a head- 
to-head configuration. In between them, a region of random 
magnetization (of length l T ) has been included. 


stant or using different initial conditions for the relax¬ 
ation process. Our results underline the crucial role of 
topological defects for physics of DWs in soft strips, and 
that of micromagnetic simulations for finding the true 
equilibrium DW structure. 


II. MICROMAGNETIC SIMULATIONS 

The system studied is a Permalloy strip of width w and 
thickness A 2 , satisfying A 2 <C re, see Fig. [l] (a). In the 
micromagnetic simulations, magnetic charges are com¬ 
pensated on the left and right ends of the strip, to mimic 
an infinitely long strip; the actual simulated length sat¬ 
isfies l > 4 w for all cases considered. The initial state 
from which the relaxation towards a stable DW struc¬ 
ture starts is an in-plane head-to-head domain struc¬ 
ture, with a region of random magnetization of length 
l r in the middle of the sample, see Fig. [l] (b). If not 
specified otherwise, we consider l r = 2 w. Material pa¬ 
rameters of Permalloy are used, i.e. saturation magne¬ 
tization M s = 860 x 10 3 A/m and exchange constant 
A ex = 13 x 10 -12 J/m. The typical Gilbert damping 
constant for Permalloy is a = 0.01, but here we analyze 
also the influence of a on the relaxation process, and thus 
consider also other values. For simplicity, we set the tem¬ 
perature T to zero, and focus on the ideal case of strips 
free of any structural disorder or impurities. 

The simulations are performed using the GPU- 
accelerated micromagnetic code MuMaxT^HUS offering 
a significant speedup as compared to CPU codes for the 
large system sizes we consider here. To calculate the mag¬ 
netization dyna mics o f the system, the Landau-Lifshitz- 
Gilbert equatioiP^H 

dm/dt = 7 H e ff x m + am x dm/dt, (1) 

is solved numerically. Here, m is the magnetization, 7 
the gyromagnetic ratio, and H e ff the effective field, with 




FIG. 2. (color online) An example of the temporal evolution 
of the relaxation, with w — 420 nm, A z — 10 nm, a = 3 and 
l T = 2 w. Relaxation towards the equilibrum DW structure 
(here, a VW) takes place via coarsening dynamics of the de¬ 
fect structure in the magnetic texture. The colorwheel in the 
middle shows the mapping between magnetization directions 
and colors. 


contributions due to exchange, Zeeman, and demagne¬ 
tizing energies. The size of the discretization cell used 
depends on the system size, but is always bounded by 
the exchange length, A = (2A//ioM 2 ) 1 / 2 « 5 nm, in the 
in-plane directions, and equals A 2 in the the out-of-plane 
direction. 


III. RESULTS 

We start by considering the effect of varying a and l r 
on the relaxation process. Fig. [2] shows an example of 
the time evolution of m(r, t) for w = 420 nm, A^ = 10 
nm, a = 3 and l r = 2 w. The initially random magne¬ 
tization evolves via coarsening of the defect structure of 
the magnetization texture towards the stable DW (here, 
a VW). During the relaxation, the total energy E of the 
system decreases in a manner that for a given geometry 
(w and A z ) depends on both a and Z r , see Fig. [3] (a) and 

(b) where w = 5120 nm and A^ = 20 nm, is considered. 
For instance, E decreases faster for an intermediate a 
[Fig. | (a)]. We attribute this behavior to the balance 
between inertial effects related to precession favored by 
a small a, helping to overcome energy barriers, and the 
higher rate of energy dissipation due to a large a. Thus, 
the relaxation time to reach a (meta)stable DW struc¬ 
ture depends on a. Fig. [3] (b) illustrates that for a fixed 
a, systems with a larger l r relax more slowly. Fig. [3] 

(c) shows that on average, the early-time relaxation of 
E towards its final value Ef exhibits temporal power-law 
decay, ( E — Ef) oc with /3 « 1.3 for the a = 0.3 case 
shown, possibly related to collective effects due to inter¬ 
actions between several topological defects during early 
stages of relaxation (Fig. [2|. 

In general, the final (meta)stable DW structure may 
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FIG. 3. (color online) The energy E(t) as a function of time 
t for w = 5120 nm and A z = 20 nm. (a) For different values 
of a and l r = 2w. (b) For different values of l r , and a = 0.3 
[resulting in the fastest relaxation in (a)], (c) shows that on 
average, the early time decay of E(t) towards its final value 
Ef obeys (E(t) — Ef) oc t _/3 . For the a = 0.3 case shown here, 
/3 ~ 1.3. Empty (filled) symbols in (c) correspond to w — 420, 
A z — 10 nm (w — 860, A 2 = 20 nm). 



FIG. 4. (color online) Phase diagram of the equilibrium 
DW structure in Permalloy strips of various thicknesses (from 
A 2 = 5 to 25 nm) and widths ranging from w — 120 nm up to 
5120 nm. The symbols correspond to observations of the vari¬ 
ous equilibrium DW structures, with phase boundaries shown 
as solid lines. Examples of the DW structures corresponding 
to the 5 different phases are shown in Fig. [5] 



depend on the realization of the random initial state. 
Thus, we consider 21 realizations of the initial random 
magnetization for each w and A z , and compare the en¬ 
ergies of the resulting relaxed configurations. The struc¬ 
ture with the lowest energy is chosen as the equilibrium 
structure, while others with higher energy are metastable 
states. Although, as discussed above, the relaxation 
times depend on a and Z r , the equilibrium DW struc¬ 
ture is found to be independent of a and l r in the range 
considered, i.e. a G [0.01,3] and l T G [w,3w\. Thus, in 
what follows, we will use a = 3 and l Y = 2w. 

The main results of this paper are summarized in 
Figs. [4] and [5j showing the phase diagram of the equi¬ 
librium DW structures for w ranging from 120 to 5120 
nm, and A z from 5 to 25 nm, and examples of these 
structures, respectively. For small re, we recover the pre¬ 
vious results^MS i. e . phases correspoding to TW, ATW 
and VW, shown in Fig. [5] (a), (b) and (c), respectively. 

For larger strip widths (w approaching or exceeding 
1/rm, depending on A 2 , see Fig.|4|, a new equilibrium mi- 
cromagnetic DW structure, a double vortex wall (DVW), 
is observed. This structure consists of two vortices with 
opposite sense of rotation of the magnetization around 
the vortex core, see Fig. [5] (d). At the phase boundary 
(blue triangle symbols pointing left in Fig. [4]), VW and 
DVW have the same energy. The DVW phase spans a 
relatively large area within the (w,A z ) space, highlight¬ 
ing the robustness of our results. 

In addition, a second new phase, with a triple-vortex 
wall (TVW) as the equilibrium structure [see Fig. [5] (e)], 
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FIG. 5. (color online) Examples of the different equilibrium 
micromagnetic DW structures: (a) TW for w = 120 nm and 
A 2 = 5 nm, (b) ATW for w — 160 nm and A 2 — 10 nm, (c) 
VW for w — 640 nm and A z — 15 nm, (d) DVW for w — 2560 
nm and A z — 20 nm, and (e) TVW for w = 5120 nm and 
A z = 25 nm. The colorwheel (top left) shows the mapping 
between magnetization directions and colors. 


is found for the very largest system sizes we have been 
able to simulate. The middle vortex of the TVW has an 
opposite sense of rotation to the other two. For w = 5120 
and A^ = 25 nm, DVW and TVW have the same energy 
(the cyan square symbol in the top right corner of Fig.[4|, 
suggesting the presence of a phase boundary between the 
two structures. Indeed, by performing a set of 10 ad¬ 
ditional simulations with w = 6144 and A^ = 25 nm 
(i.e. outside the phase diagram in Fig. |4|, suggests that 
TVW is the equilibrium DW structure for very large strip 
widths. This structure has been observed in experiments 
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FIG. 6. (color online) Examples of metastable DW structures 
observed for a system with w — 5120 nm and different thick¬ 
nesses A z : (a) Two vortices and an antivortex (2V+AV), 
= 25 nm, (b) and (c) three vortices and an antivortex 
(3V+AV), A z = 5 nm, (d) three vortices and two antivortices 
(3V+2AV), = 10 nm, (e) four vortices and an antivortex 

(4V+AV), A^ = 5 nm, (g) four vortices and three antivor¬ 
tices (4V+3AV), A z = 5 nm, and (h) five vortices and two 
antivortices (5V+2AV), A z = 5 nm. 


as a metastable state for smaller systems 28 Notice also 
that the middle part of the TVW [Fig. [5] (e)], exhibiting 
four line-like 90° DWs meeting at a vortex core in the 
middle of the TVW, resembles the typical Landau flux- 
closure magnetization patterns observed for rectangular 
Permalloy thin films 3 ^®. 

Following the relaxation from a random magnetization 
initial state, the system may in general end up into var¬ 
ious metastable states with higher energy than that of 
the equilibrium DW. Sometimes these metastable states 
have even a higher probability than the equilibrium one, 
estimated here from the sample of 21 relaxed configura¬ 
tions. Fig. [6] shows some of the metastable states found 
for a large strip with w = 5120 nm and different values of 
A 2 ; for strips with smaller lateral dimensions, different 
metastable states tend to be less numerous and have a 
simpler structure. Despite their apparent complexity, all 
the metastable DW structures shown in Fig. [6] respect the 
basic principles of topology of DWs. Each of the DWs 
are composed of topological defects, with an associated 
winding number: +1 for vortices, -1 for antivortices, and 
±1/2 for edge defectIn a DW all the topological 
defects have to be compensated, i.e. the total winding 
number is equal to zero. In the case of the DVW, the 
two topological vortex defects are compensated by four 
edge defects [Fig. [5] (d)]. For the metastable state of two 


vortices (with the same sense of rotation) and an antivor¬ 
tex [2V±AV, see Fig. [6] (a)], there are two vortices and 
only two edge defects. Thus, in order to compensate the 
topological defects, also an antivortex appears. In gen¬ 
eral, we have observed that in a DW with N vortices with 
the same sense of rotation, there must be N — 1 antivor¬ 
tices to get a zero total winding number, see Fig. [6] (a) 
and (d) for examples with 2V±AV and 3V±2AV config¬ 
urations, respectively. When some of the vortices have 
oppposite sense of rotation, more complex scenarios are 
encountered, with examples shown in Figs. [6] (b), (c), (e), 
(f), (g) and (h). Notice also that two DW structures with 
the same elements can look very different, see e.g. the 
two 3V±AV DWs shown in Figs. [6] (b) and (c). All the 
DW structures found, both the equilibrium ones in Fig. [5] 
and the metastable states in Fig. [6] obey the principle of 
compensation of topological defects to yield a total wind¬ 
ing number of zero. The richness of the equilibrium phase 
diagram and the large collection of metastable states in¬ 
dicate that for wide/thick strips in particular, the micro- 
magnetic energy landscape is quite complex, with a large 
number of local minima. This is also in agreement with 
our observations of power-law energy relaxation. 


IV. SUMMARY AND CONCLUSIONS 

To summarize, we have performed an extensive set 
of micromagnetic simulations to study the equilibrium 
and metastable DW structures in Permalloy strips of a 
wide range of widths and thicknesses, as well as the re¬ 
laxation dynamics starting from random magnetization 
initial states. The general trend of our results is that 
both the equilibrium and metastable DW configurations 
become increasingly complex (i.e. they consist of an in¬ 
creasing number of topological defects) as the lateral strip 
dimensions increase. We note that somewhat analogous 
behaviour - i.e. existence of equilibrium magnetization 
configurations with increasing complexity as the system 
size increases - is observed also in some other systems 
such as three-dimensional cylindrical elements with per¬ 
pendicular anisotrop}®^. 

Several remarks are in order: first, for strips with 
even larger lateral dimensions one may in principle ex¬ 
pect more complex DW patterns - possibly with four or 
more vortices with alternating sense of rotation. These, 
however, are currently beyond the reach of our available 
computing resources. Second, our phase diagram allows 
one to check if experimental observations of the various 
DW structures in wider strips are equilibrium configura¬ 
tions or metastable states. According to our review of 
the experimental literature, most observati ons o f DVWs 
and TVWs appear to be metastable stateeP^^. Third, 
while the equilibrium structures we find are certainly sta¬ 
ble in the absence of external perturbations such as ap¬ 
plied magnetic fields, it remains to be seen how their field 
driven dynamics is like, and whether wide strips with a 
relatively weak shape ansitropy are able to support the 
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DWs as compact objects also when external perturba- ACKNOWLEDGMENTS 
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